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Abstract 

We study renormalization of Coulomb-gauge QCD within the Lagrangian, second- 
order, formalism. We derive a Ward identity and the Zinn- Justin equation, and, with 

_ , the help of the latter, we give a proof of algebraic renormalizability of the theory. 

^ ■ Through diagrammatic analysis, we show that, in the strict Coulomb gauge, g'^D'^^ 



O 

o 



^ 



is invariant under renormalization. (D is the time-time component of the gluon 



— ' 

^^ . propagator.) 



1 Introduction 

In nonabelian gauge theories, among a variety of gauge choices, the Coulomb gauge 
is one of the most important ones. The theories with this gauge are described in 
terms of physical fields, so that the unitarity is manifest. Within the Hamiltonian, 
first-order, formalism, formal or algebraic renormalizability of the Coulomb gauge 
has been studied in [1]. Since then, introducing an interpolating gauge, which inter- 
polate between a covariant gauge and the Coulomb gauge, Baulieu and Zwanziger 
have proved algebraic renormalizability of the theory [2]. In taking the Coulomb- 
gauge limit, a phase-spaces representation is used there. Despite its importance, the 
proof of renormalizability of the Coulomb gauge within the Lagrangian, second- order, 
formalism per se is still lacking, to which the present paper is devoted. 

In this paper, we are concerned about Coulomb gauge QCD, whose Lagrangian 
density is obtained by adding the gauge fixing term, £g.f. = ^{diAD{djAi)/0^, to 
the QCD Lagrangian density. In the Coulomb gauge, there is an inherent problem 
of appearance of 'energy divergences', which are characteristic of an instantaneous 
closed ghost- and Ag-loops: 

fpF(p....), (1) 

J Zn 

where F is independent oipo, the temporal component of the four- vector P^ = [po, p) 
and "..." indicates a set of external momenta. There also appear ill-defined integrals 
of the forms, 

Glp, ...), (2) 



2vr pl-p^ + iO+ 

dpQ dqo Po go 



iJ(p,q,...), (3) 



where G (H) is independent of po {po and go) • 

A numerous work has been devoted to the energy divergence problem and, by 
now, the following results are established. 

1) In the Hamiltonian, phase-space, first order, formalism, energy divergences like 

(1) do not appear [3] in the first place. 

2) Using the correspondence formula which equates amplitudes in a covariant gauge 

to those in a gauge without ghosts, Cheng and Tsai [4] "indirectly" showed 



that, when all relevant contributions are added, cancellation occurs between 
the energy divergences (1). 

3) With the help of an interpolating gauge, which interpolates between a covariant 

gauge and the Coulomb gauge, in the phase-space formalism, it has been shown 
by Baulieu and Zwanziger [2] that the cancellation occurs between different 
contributions which turn out to energy-divergent ones in the Coulomb-gauge 
limit (see, also, [3]). 

4) Ill-defined integrals like (2) can be set equal to zero [5] and another type of ill- 

defined integrals, Eq. (3), are connected [3, 5] with the so-called Vi + V2 terms of 
Christ and Lee [6], which arise through correct treatment of operator ordering 
in the Hamiltonian. 

5) It has been shown in [7] that the cancellation of energy divergences and renormal- 

izability is compatible in an example in which quark-loop subgraphs are inserted 
into the second-order gluon self-energy graphs: As mentioned in 1) above, in 
the phase-space formalism, two integrals over the internal energies converge. 
However, in relation to renormalization, energy-divergences re-appear. Thanks 
to the Ward identity, these energy-divergent contributions cancel out. 

Therefore, when perturbative computations in the present Lagrangian formalism are 
properly handled, cancellation should occur between the energy divergences. We are 
not concerned, in this paper, with the energy- divergence issue anymore. 

As stated above, we are are interested in the Coulomb gauge QCD. Nevertheless, 
we proceed, as far as possible, with more general gauge choice, £g.f. = ^{.diA\Y / {2q.) 
with arbitrary a, which is usually called the Coulomb gauge. The "genuine" Coulomb 
gauge, which is obtained by taking the limit a ^ 0^, is called the strict Coulomb 
gauge. Incidentally, in the case of a 7^ Coulomb gauge, when compared to the strict 
Coulomb gauge, much worse energy divergences arise (cf. Eq. (5) below), so that, at 
the present stage, such a gauge is useless for practical perturbative calculations. 

In §2, we present a Ward identity and the Zinn- Justin equation, which are derived 
in Appendix A. In §3, on the basis of the Zinn- Justin equation, we construct renormal- 
ization counterterms in a recursive way. In §4, we prove an algebraic renormalizability 
of the theory and obtain the identities between the renormalization constants. With 



the aid of the Ward identity, we obtain an additional identity Z\ = Z12/Z31 (Eq. 
(41)). We then show that, in the strict Coulomb gauge, g'^D^^ {D^^ is the time-time 
component of the gluon propagator) is invariant under renormalization. In [1, 2], this 
proposition is proved using the Ward identity. In contrast, in the present Lagrangian 
formalism, use of the Zinn- Justin equation is sufficient for proving this. For illus- 
trative purpose, the forms of some one-loop renormalization constants in the strict 
Coulomb gauge are also displayed. §5 is devoted to summary. Appendix A briefly de- 
scribes a derivation of the Ward identity and the Zinn- Justin equation. In Appendix 
B, using the identity that is dereived from the Zinn- Justin equation, we obtain, in 
the strict Coulomb gauge, an additional identities Zi = D' = 1, Eq. (B.5). In Ap- 
pendix C, we briefly describe a derivation of some one-loop renormalization constants 
in the strict Coulomb gauge. Some formal diagrammatic analyses are summarized in 
Appendix D. In particular, we show that, in the strict Coulomb gauge, r^^T^-vertex 
(cf. Eq. (4)) is not a renormalization part, so that Zi = 1 in the minimal subtraction 
scheme in dimensional regularization. 

2 Ward identity and Zinn-Justin equation 

As the content of this section is standard, we describe briefly. Greek indices fi, u, ... 
run over 0, 1, 2, 3, while Latin indices i, j run over 1, 2, 3. We use P'^ for denoting a 
four vector P^^ = {p^, p) and p' for denoting a three vector. 

The effective Lagrangian density of Coulomb-gauge QCD with one quark flavor 
(generalization to the case of several quarks is straightforward) reads 

£cff = Ccs-—J^a^a [^:Fa[A; x] = diAKx)^ , 

Cos = --FrFaf,, + ^{t0 -m-gtaAa)^ + dir]aDl,{x)7]b, (4) 

where F^ = 9^^^ - 9M;^ - gfabcA^b^'',, ta = A„/2, and Z^^,(x) = dabd^' + gfabcA^d^)- 
Generalization to other nonabelian gauge theories is straightforward. In the La- 
grangian formalism adopted here, the fields propagators can be extracted from the 
bilinear (with respect to the fields) terms of Ccs in Eq- (4). For the purpose of 
later use, among the propagators, we only display the forms of the gluon propagator 
A'^'^{Q) (see, also, [8, 9]) and FP-ghost propagator A{Q), together with f]a{Q)Alric- 



vertex factor V^bciQ)'- 






A'^'^(g) = g',g^, ^, ..? +-^-^^^-1^, (5) 



A(Q) = -^, (6) 

VL(Q) = ^/aftcgV (7) 

Propagators are diagonal in color space, so that the color indices are suppressed. 

The gauge-field part and the quark part of £eflf is invariant under the infinitesimal 
gauge transformation: A^ -^ A'^+D'^^{x)eij{x), ip -^ ip—igeataip, and-?/' -^ ip+igipeata- 
Ces is invariant under the BRST transformation of the fields 

6A^,{x) = D,^,(xHC = sA^C, (8) 

SVaix) = --fabcVbix)r]c{x)C = sr]aC, (9) 

Stp = igriatai^C = #C , ^'4^ = igi>VataC = si'C, (10) 

Sf]a{x) = -J^a[A;x]( = Sf]aC. (11) 

a 
Here ( is an x-independent infinitesimal Grassmann number with the same ghost 
number as fja- C^s is not invariant under the Lorentz transformation but is invariant 
under the spatial rotation. Then, we treat the spatial component A^ and the temporal 
component A^ of A'^ separately. We let x" run over the fields, A^, A°, r]a, ip, and 
■ijj, and introduce a set of c-number external sources Kn{x) that couple to sx"[X5 2;], 

Quantum effective action F is defined by the following implicit functional integro- 
differential equation: 



{x--r)--^iv'.-v.) 



(12) 



5x" ^ ^ 5r]a 

where '/?' denotes right differentiation. We show in Appendix A that 

r[x, fj, K] depends on fj only through K^ + d'fja ■ (13) 

We introduce F through 

f = T + ^jdS:F'[A-y]. (14) 



The Ward identity, which is derived in Appendix A, reads 



(fx 



6rT 6lT 



{x) = do d^x 



Va^-^ + K^.^'^ 



MO 



5K. 



Va 



X 



(15) 



where 'L' denotes left differentiation. Integration over xq yields the Zinn- Justin equa- 
tion. 







d X 



SrT 6lT 



X] 



The leading term (in the loop expansion) of F is 



r, r . 



ro[x, V,K]= / d'^x (£cflf[x, ^, ; a;] + i^„sx"[x; x\) , 



(16) 



(17) 



which is invariant under the transformation (8) - (10). It should be noted that the 
fields x" and fj in Fq are the renormalized ones. 



3 Recursive construction of counterterms 

In this section, we construct the renormalization counterterms that preserve the sym- 
metry condition (16). We follow the procedure in standard text books [10, 11], so 
that we briefly describe. 

Let us use dimensional regularization by continuing spacetime dimensions from 
4 to d. We employ the loop expansion for F, F = Y1,'n=q'^n and decompose Fat as 



N 



^(oo) 



^(oo) 



^RN + Tat , where Fjy is the ultra-violet (UV) divergent contribution, i.e., it 



diverges in the limit d —^ A. We adopt the minimal subtraction scheme. 



The symmetry condition (16) leads to Y,n'=q {J^n', '^n-n') = (A^ = 0, 1, 2,...). 
We proceed in a recursive way. We assume that, for all M < A^ — 1, all UV-divergent 
contributions from M-loop diagrams have been cancelled by counterterms T m {M 
= 1, 2, ..., A^ — 1). (This recursive procedure is justified a posteriori in §4.) Then UV 
infinities can appear only in the A^' = and N' = N terms, and the infinite part of 

= 0. Here, Fq is as in Eq. (17). 



the above condition reads (Fn, F 
We write T}^' as 



,(oo) 

N 



+ r 



^(oo) 
N 1 



^^n\x.V.K] 



d^x (In[x, a + KnVlilx; x 



where In and V^ are local functions of x" and rj and their derivatives [11]. For the 
time being, we drop the suffix A^. Substitution of Eqs. (17) and (18) gives 



,4 I dRCcs[x,V;x]^n! 1 , dRl[x,v;x] 



'^ [ a^q;; ^'fe ^1 + ^#r'''^"i« ^1 j - ° • <'^) 



,4 /''9ijsx"[x;a;]^„. _ 1 , dRV"'[x;x] 



From here on, we follow the procedure in Sec. 17 of [11]. We introduce 

f ^'^ [X, ^] = J ^^^ {^os[x, V; x] + el[x, v; xfj 

^ Jd^xC^^^ix^mx], (21) 

A«"[x;x] = sr'[x;x]+eV-[x;x], (22) 

with e infinitesimal for technical reason. Then, Eqs. (19) says that T^'^' is invariant 
under the transformation, 

x"(x)^x"(:^) + A(^)"[x;x]C, (23) 

while Eq. (20) tells us that this transformation is nilpotent. 
The most general form of the transformation (23) is 



Al, -^ A\ + 



a 



A° ^ A 



a 



B;S'^\-gD^£Alr^^C, 



9 77(f) /■ 
ip -^ ip + igriaT^'^-^C , 

where bI^^\ BJ^^ , D^l, D^^l, and E^l (= -^i'i) are constants, and T^'^ is some 
matrix acting on the quark field. Here D^^jj^ = iD^abc)N = fate + e(£'a6c)Af, T^"-* = 
{Tlf^)iy = ta + eiTa)N, etc. Even if we introduce different T^^^^'s for different compo- 
nents of ip, we will have the same result as the one obtained below. Imposing the 
condition of nilpotence, we obtain E^^l = D^^l = D^^ = C^'^fabc, B^^ = {CD)^^6ab, 
B^J = {CD'Y'Uab, and T^^) = C^'Ha, where C^'\ {CDY'\ and {CD'Y'^ are some 
constants: C^'^ = C^,f = 1 + eCjv, (CD)^'^ = {CD)%^ = 1 + eEJ^'=oCiv'^iv-7V' 
(Co = Dq = 1), etc. {Cm and Dm {M < N) have already been determined at 
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lower-order stages.) Thus, we have, with obvious notation. 



K - 


- a: + 


CD {d% - gfabcAlr],) C 


1 


K - 


- K + 


CD {{D'/D)d%a - gfabcAtvc} C 


Va - 


^ Va + 


CD (--gfabcVbVcC) 




ij - 


-^ ij + [CD{tg7]J^,pC)f'^ , 







(^) 



(24) 
(25) 

(26) 

(27) 



where g = g/D. Introducing A° = {D/D') A^, one can rewrite Eq. (25) as 



i° ^ i° + \CD (d% - -gfabcA^Vc) C 



(^) 



C^^^ in Eq. (21) must be invariant under all the symmetries of the original La- 
grangian, i.e., (1) invariance under the spatial rotation, (2) global gauge invariance, 
(3) antighost translation invariance, (4) ghost number conservation, and (5) quark 
number conservation. Furthermore, as stated in (13), F does not involve dofj. (For 
completeness, this is proved in Appendix D through formal diagrammatic analyses.) 
Taking all of these facts, we find that the most general renormalizable interaction 
takes the form: 



^(e) 



£cff + e/ 



4a + ZJ^^\d,naWVa) - gdl:l{d,na)AlVc , 



(2^ 



where Z^'^\ d^^^ are unknown constants, and Cl\ is the renormalizable term that 
involves only the quark and gauge fields. 

Imposition of the invariance under the transformation (24) - (27) yields 



£(.) 






(^) 



(29) 



The tilde on F'^J indicates that the field strength is to be calculated using g for g^ 
while the tilde on F^* indicates that the field strength is to be calculated using g for 
g and A'^ for A*^. All the constants appearing above are of the form, Z^\ = (Z^l)]^ = 
1 + e(Z3i)^, (Zsi/DY^^ = (Zai/D)!^^ = 1 + eE%'=o{Z3i)N'D],lj,, {D,' = 1), etc. 



For obtaining further informations on the constants appearing in Eq. (29), we 
make following observations: 

(i) Gauge-fixing independence of the locations of the physical poles (zeros) in the 
transverse-gluon- and the quark-propagators (two-point effective actions). 
Kobes, Kunstattar, and Rebhan [12] showed this proposition to hold by using 
a set of identities that determine the gauge dependence of the effective action. 

(ii) Argument in [12] goes as it is even in rf-dimensional spacetime. 

(iii) In the case of covariant gauge, the poles of the transverse-gluon propagator 
are at po = ±p (P^ = 0) and the poles of the quark propagator are at po = 
±(p^ -|- ^pjj)^/^ (P^ = rrip^), where mph is the physical quark mass. Thus, the 
dispersion relation of the transverse-gluon (quark) mode is Lorentz invariant. 

Our Z^'^^^s are the ones that are computed in the minimal subtraction scheme in 
dimensional regularization. For taking the facts (i) - (iii) into account, the on-shell 
renormalization scheme (RS) should be employed. Let (^3^^)011, etc. be the on-shell 
RS counterpart of Z^\ , etc. Then, (i) - (iii) and the on-shell RS counterparts of Eq. 
(29) tell us that {Z^^U = (^if )on and {Z^fU = (4?)on. Since (zif),,, - Z^ , etc. 
are UV finite, we obtain^ 

yi'^) _ y(<^Y y(<^) _ yi'^) (— yi'^)\ ('3r\\ 



4 Algebraic renormalizability and the form of r*^°^) 

In this section, we show the algebraic renormalizability of the Coulomb-gauge QCD. 
On the basis of the recursive construction in the last section, we deduce the form of 
V'^°°\ Reviving the suffix A^, setting e = 1, and summing over iV, we obtain from Eqs. 
(28) and (29) with Eq. (30), 

00 

N=l 



^For arriving at the strict Coulomb gauge (a ^ 0), one might encounter a subtle limiting proce- 
dure. The proposition (i) has been proved within the interpolating gauge in [2]. 
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~2 



'' ^d^Aa.d^Ai - d^Aid.Al) + Z,2doAld^A'i - -H5,.a09M° 



+9ZufabMA])AlAi + (7Zi2/a6c5°A^,A°A„ + (^Zia/afecSM^ A,,A° 

Q Q 

~r ZilJ ahcj adeAuAcjAiji^A^ — Z 42 J abc Jade A)j Ad A^A^ 

+Zs{dif]a)d'r]a - gZifabc{dif]a)Alr]c 

+Z2i> {i0 - (m - 5m)) ^ - gZ^ni'taAi-fj^p - gZ^ui'taA'^^-f'^^p , (31) 

where, with obvious notation, Z^i = 1 + Z]^=i(^3i)Af, etc., and 

Z32 = DZsi/D', Z33 = D^Zsi/D'^ Zn = Z,,/D, Z^2 = Z,^/D', 

Zi3 = DZsi/D'\ Z41 = Zsi/D\ Z42 = Z,i/D'^ Z, = Z^/D, 

Z^ii = Z2/D, Z^i2 = Z2/D', 5m = (l- Z'„jZ2)m . 

Here, D = 1 + J2n=i ^n, etc. From these relations we obtain the identities: 

Z31 _ -^3 _ ^11 _ Z32 _ ^33 _ Z2 ,_ . 

Z\\ Zx Z41 Z12 Zi3 Z^ii 

= D') . (32) 



^31 Z12 Z32 Z2 



AJ _ 7^/^AJ 
^Ba — ^31 ^a ) 


^Ba 


IpB = Zl^^lp , 


i>B -- 


m R = m — 6m . 





Z\2 Z42 Zi3 Z^i2 

Recalhng Eq. (14), we finally have, for the density of (Fq + F*^°°M , 

i = i - ^{d.A:){d,Ai) , (33) 

where a = Z^ia. We now define the bare fields and parameters according to 

_ 7I/240 ^ _ 71/2^ ^ _ 7I/2- 

- ^33 ^a > VB - ^3 V^ VB - ^3 V, 

Zy^i) , gB = ZnZ^x^^'^g , Ob = ^31" , (34) 



Then, from Eqs. (31) - (33) we see that i may be regarded [10] as the initial La- 
grangian density Ccs{AB,riB,f]B,4'B,'4'B'igB,iTT'B,C(B) written in terms of bare quan- 
tities. 

By construction, the density i in Eq. (31) is invariant under the transformation 
X" ^ X" + ^X" that is obtained from Eqs. (24) - (27) by setting e = 1 and summing 
over N, which reads, in obvious notation, 

AAl = SAi^C = C (Dd'v, - gfatcAlVc) C = ^Z^^' (d'VBa - gBfabcA^BtVBc) C 
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C Z\\Z. 



1/2 



12 



d^riBa - gBhhcAs^riBc] C 



Ar/„ 



Aijj 



ZllsA^BaC , 



SriaC = C ( -^gfabcVbVc ] C = ^^sf 



c 
zl 



-^gBfabcVBbVBc ) C = Zr^^STJBaC , 



S^C = C (tgVata^) C = |- (^^ 



1/2 



(igVBatai^B) C = Z^S^psC , (35) 



where use has been made oi g = g/D and Eq. (34). Then, from Eq. (18) with Eqs. 
(22), (23), and (35) we have 



f + E r^^ 



N 



N=l 



d X 



+ KnSx'' 



(36) 



d'^X (£cff[XB, Vb] + ^BnSXBlXB^ 



where we have used £ = CcsIxbiVb] (cf- after Eq. (34)). On the right-hand side of 
Eq. (36), K* and d'^'f]a appear in the form. 



CZ^'Kl + {d'f^a)\ \^^'na ' ~gfabcAuVc] ■ 

Recalhng here the property (13), we obtain the relation C = Zi. Thus, 

{Z^, Z^, Zn^, Z^,Zx) cannot be arbitrarily chosen but are determined uniquely: 



(37) 



Z 

z. 



Va 



a/2 



7I/2 

"31 ! 



^a — ^11^3 /^12 , 



^31^3/^2 



1/2 



(3^ 



1/2 



From these relations and Eq. (34), we find KbtiXb ~ (^31^3) -^nX"? where sum- 
mation over n is not taken (n = 1, 2, .., 5). Note that the factor (Z31Z3 



IS common 

for all n. This fact, together with the relation i = Ccs{Ab, tjb, tjBi "^b, "^b] Qb, itT'B, <^b) 
obtained above, justifies [10] a posteriori the recursive procedure in §3. 

It should be emphasized that we have uniquely deduced the results (38), which is 
in contrast with those in [2]. 



Derivation of Zi = ^12/^31 

Following [2] , we start with the Ward identity written in terms of the bare quantities: 
5rTb SlTb 



d-^x 



6xb ^Kbu 



[X 



do / d X 



SrTb , r^ h^B 



"SA%^ 



Ba 5K, 



VBa. 



\X 



(39) 
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According to our results on renormalization, the effective action Tb is finite [2] when 
expressed in terms of renormahzed quantities, 

f b(xb, Vb, Kb; Qb, rriB, Ob) = f (x, V, K; g, m, a) . 
Through this change of quantities, Eq. (39) turns out to 



d'^x 



(5x" 5Kr, 



I N ^31^1 

[x] 



do I (fx 






5AI ''-SK,^, 



{x) . (40) 



^12 

Since the left-hand side of this equation is UV finite, Z^iZi/Zi2 must be finite. This 
implies that in the recursive procedure described in §3, the UV-divergent part of 
Z^iZi is equal to that of Z12 in each loop order. Then, we obtain the relation 

where use has been made of Eq. (32). 

Strict Coulomb gauge 

In Appendices B and D, we show that, in the strict Coulomb gauge, Zi = D' = 1, 
which is in accord with Eq. (41). Using this in Eq. (32), we obtain 

Z12 = Zn2 = ^31 ? ^13 = ^32 ) ^vi2 — ^2 • 
Then, the relation qbA^^ = gA^ (cf. after Eq. (34)) and then 

glDl' = g'D'^ , (42) 

hold, where D^^ is the time-time component of the gluon propagator. 

It is worth mentioning here that, in the Hamiltonian, first-order, formalism [2], 
the identity (42) is deduced from the Ward identity. In contrast, as seen above, for 
deriving Eq. (42) in the Lagrangian formalism dealt with here, use of the Ward 
identity is not necessary. As a matter of fact, D' = 1 and Zi = 1 plays a role. 
The former is derived from the Zinn- Justin equation in Appendix B and the latter 
is obtained in Appendix D through formal diagrammatic analysis as in the case of 
Landau gauge. Physical importance of the relation (42) is fully discussed in [1, 2]. 

As a matter of course the Ward identity contains richer informations than the 
Zinn- Justin equation, so that different relations are to be obtained between various 
(renormahzed) amplitudes. 
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One-loop wave-function renormalization constants and dm in the strict 
Coulomb gauge 

As an illustration, we display here the results for the wave-function renormalization 
constants and 6m to one-loop order in the strict Coulomb gauge (see Appendix C): 

_ 7^2 1 9^ I 

'^ - &dh^ (^^) 

9^ 1 

ivr^ a — 4 

The second terms on the right-hand sides of Eqs. (43) - (45) come from the quark-loop 
diagram. These Z's and Z3 satisfy the identities (32). 

We explicitly confirmed in the Lagrangian formalism, as in other literatures in 
different formalisms, the absence of energy divergence in one-loop order (Appendix 

C). 

5 Summary and discussion 

In this paper we have addressed the problem of renormalizability of the Coulomb- 
gauge QCD within the Lagrangian, second-order, formalism. Starting with the Zinn- 
Justin equation and following the procedure as in [11], we have proved a formal or 
algebraic renormalizability. The renormalization constants for the external sources 
-ft'*, K^, K^^, K^, and K^, which couple, respectively, to the composite operators 
sAl^, sA°, srja, sip, and sip are uniquely determined. 

We have derived the Ward identity (15) and from which we have obtained the 
identity Z3 = Z32/Z31. With the aid of the Zinn- Justin equation, we have shown that 
gA^ is unchanged under renormalization. In [1, 2], this is derived using the Ward 
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identity. In contrast, in the present Lagrangian formalism, use of the Zinn- Justin 
equation is sufficient for proving this. 
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A Ward identity and Zinn-Justin equation 

fja-dependence of F 

Changing the integration variable fj'^ ^ fi'^ + Sfj'g^ in Eq. (12), we obtain 



= Jv{x',v')6v'a 
'Vix',v')6r]', 
From this equation, we get 



ADiMH ' ^""^ 



e 



/A(...) 



i6Ki Sr]a 






This relation tells us that F depends on fj only through Ki + d^f]a. 

Derivations of Ward identity and Zinn-Justin equation 

We introduce infinitesimal variations [1, 2] 

X'"(:^) - x'' + f{xo)sx'\x)C (n = l,...,5), (A.l) 

n'Ax) ^ n'a{^) + f{x,)st{x)C, (A.2) 

where sx " and sf]'^ are as in Eqs. (8) - (11), and /(xq) is a xo-dependent function. 
When / is a constant, it reduces to the BRST transformation. Making this change 
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of variables in Eq. (12), we obtain 



lT^ix',v') 



X exp 






5^a 



ijl'a - Va) 



0, 



where p is the BRST charge: 



(A.3) 



P[X] = -Fj^'sA',, + z^p^sij' + K^sr^', 
We rewrite the term being proportional to pd^f as 



v{^,n' 



5rT 



ijd*y(...) 



i J d^y {^Cesix', v') + KnSX 



X exp 



n7^2 "a. 






(A.4) 



Carrying out the partial integration, we obtain 



E,. (A.4) ^ - jvi^.f^ [j,^. (,:M + ;,;0^) ao/ 



,i/d4y{...) 



(A.5) 



From the definition of rja, we can make a replacement 1]'^ -^ f]a ■ 

From here, we proceed as follows: i) Substitute Eq. (A.5) in Eq. (A.3), ii) 
translate to F by Eq. (14), and iii) integrate by part with respect to Xq for the term 
that involves dof. Then, using the arbitrariness of /{xq), we finally obtain the Ward 
identity (15), 



d-^x 



5rT 5lT 
5x" SKn 



(x) = do d^x 



Va 



SrT 



+ K, 



SrT 



'^L 



6 A?. "SK 



na] 



[x). 



(A.6) 



Integration of Eq. (A.6) over xq yields the Zinn- Justin equation (16). Eq. (A.6) holds 
for arbitrary a (gauge parameter). 
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B Derivation of Z31/Z12 = i^' = 1 in the strict 
Coulomb gauge 

By functionally differentiating the Zinn- Justin equation (16) with respect to A'^, A] 
and r]c, we obtain [13], after Fourier transformation, 



P^W>^{P) = W>^{P)%{P), (B.l) 

which is diagonal in color space, so that the color index is dropped. 11'^'^ is the two- 
point gluon effective action, from which the gauge-fixing term is dropped. Il^{P) is 
defined by 

^9fabM'i{x)Vc{x)Uy))tr ^ Sadtl^iP) , (B.2) 

where the suffix "tr" stands for truncation and "F.T." means to take Fourier transfor- 
mation. The FP-ghost self-energy part n(-P) is related to Il^{P) through p*nj(P) = 

n(P). 

Let us introduce a tensor decomposition of U"^: 

U^'^iP) = g^g^^g'^A + j^i^B + po(«>^ + P%'')C + n'^n^'E , (B.3) 

where n^ = (1, 0) and p- = p^ —pon^. Substitution of 11^'^, which is written in terms 
of bare quantities, into the "bare counterpart" of Eq. (B.l) yields two equations, one 
of which reads 

(p'' + IIb)Cb- [i-^^\eb = Q. (B.4) 

Another equation leads to the same result as the one obtained below. 

When expressed in terms of renormalized quantities, this equation should become 
[2] an UV-divergence free equation. Relations between the bare- and renormalized- 
quantities are Cb = (-^31-^33)" C, Eb = Z^^E, and p'^ + H^ = Z^^ \r? + 11). In 
Appendix D, we show in the strict Coulomb gauge that XIq is UV finite and Zi = 1. 
Then, 

IlBo = j-^I^ZsZU'^-^jUo = -Uo, 

where gB is gB on the left-hand side of the "bare counterpart" of Eq. (B.2). Then, 
Eq (B.4) turns out to 



(p^ + n)c^z.z,'(i-^)^. 
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Since the left-hand side is UV finite, same reasoning as in §4 (cf. Eq. (40)) leads to 

Z^=D' = 1. (B.5) 

These relations holds for the strict Coulomb gauge and are in accord with Eq. (41). 

C Derivation of Eqs. (43) - (45) 

The strict Coulomb gauge gluon propagator is given by Eq. (5) with a = 0. Straight- 
forward computation using Eqs. (5) with a = 0, (7), and the forms of 3- and 4-gluon 
vertices yields, for the one-loop UV-divergent contributions to Il^'^{Q) (Eq. (B.3)) 
(see, also, [1, 9]), 



.2 



A 



_J^ l_ 

1207r2 d-A 



15g2 + 96g2 ) - (60g2 + 96^^) + lOQ^ 



q^ 1 
" = ll^I3ll93-48-101. 

Here, each term on the right-hand side of each equation is the contribution from the 
following one-loop diagrams: 

• First term ^- The diagram that includes two transverse-gluon propagators. 

• Second term ^- The diagram that includes one transverse-gluon propagator and 
one Aq propagator and the tadpole diagram. 

• Third term ^- The diagram with a quark loop. 

From Eqs. (C.l) and (31), we extract Z's, as given in Eqs. (43) - (45) in the 
text. We briefiy describe how Eq. (C.l) is derived. The (l)-type energy-diverging 
integrations appear. As mentioned in §5, it is known that, when all relevant contribu- 
tions are added, cancellation occurs between them. Here, we substantiate this for the 
one- loop gluon two-point function in the Lagrangian formalism. Energy divergences 
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arise only in the transverse-gluon two-point function. Explicit computation of the 
energy-divergent contributions of the type (1) yields 



(n^,)^ ^ SgXtJ 



d^P 1 /^i_. pV 



5'^ 



(27r)4(k-p)2 \^ p2 y ' 

Here ~ indicates the energy-divergent contribution, and Uj {j = 1, ...,4) is the con- 
tribution from the following one-loop diagrams: 

• Hi <— The diagram that includes one transverse-gluon propagator and one Aq 
propagator. 

• 112 ^ The diagram that includes one Aq propagator (tadpole diagram). 

• lis ^ The diagram that includes two Aq propagators. 

• 114 ^ The diagram with FP- ghost loop. 

It can readily be seen that the cancellation occurs between the four contributions in 
Eq. (C.2). 

The ill-defined integrals like (2) also appear. As mentioned in §5, such integrals 
can be set equal to zero. 

D Diagrammatic analyses 

In this Appendix, we formally carry out some diagrammatic analyses. We start with 
following observations. Consider a diagram G that includes FP-ghost external lines. 

a) From Eq. (7), we see that the vertex factor for the external fjA^rj-veTtex from 
which the outgoing ghost goes out is independent of the internal loop momenta. 



b) Strict Coulomb-gauge case: From Eq. (7), the vertex factor for the external fj{P~ 
Q)A^ {Q)ri{P)-Yertex, into which the incoming ghost enters, is proportional to 
Pj — Qj. A^{Q) constitutes the gluon propagator (5), the transverse part of which 
includes 6^^ — q^q^/q'^. Then we have {pj — qj){5^^ — q^q^/q^) = Pj{S^^ — q^q'^/q'^), so 
that this vertex factor turns out to be independent of the internal loop momenta. 

UV divergent contribution to the ghost propagator n(P) 

Consider a diagram G that contributes to n(P). The vertex factor (Eq. (7)) for the 
external fi{P)A'^ri-veTtex from which the outgoing ghost goes out is proportional to 
p*. The vertex factor for another external fj{Q)A^r]-vertex, into which the incoming 
ghost enters, is proportional to q^ . Because of the rotation symmetry, we have, after 
loop-integration, n(P) oc p'^J-', where JF is a dimensionless function of Pq, p, and m. 
Since we have assumed that the divergent part of 11 is a local function of the fields 
(cf. after Eq. (18)), we have Il^°°\P) oc p'^/{d - 4)" {n = 1, 2, ...), so that £(") in Eq. 
(21) does not involve the term {dofj){dorj). 

UV finiteness of fjA^r] three-point functions 

Let J-'{P,Q) be a fjA^{Q)ri{P) three-point function. From its Lorentz structure, JF is 
of the form JF = J-'ipo + J-'2qo- This, together with the above observation a), shows 
that the degree of UV divergence is —1, so that JF is UV finite. Then, C^'^^ in Eq. 
(21) does not involve the term {dof])A^ri. 

UV finiteness of no(^) in strict Coulomb gauge 

Consider a diagram G that contributes to IIq{P), Eq. (B.2) with /i = 0. From the 
Lorentz structure no(-P) oc po- This, together with the above observations a) and b), 
shows the degree of UV divergence of the diagram is —1. 

Zi = 1 in strict Coulomb gauge 

Consider an fjA^r] three-point function in the strict Coulomb gauge. Above observa- 
tions a) and b) tell us that the fjA^r] three-point function is UV finite and then, in 
the minimal subtraction scheme, Zi = 1, as in the Landau gauge case. 
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